% just copy of FN.m
% latest date 23.10.2002

function[z] = Kill50(hz,notch),
j = sqrt(-1);

aaa = length(notch);

RQ = 1;
b = 0;
for i = 1:aaa,
b = b + sig(hz,notch(i));
end;
b = b * RQ;
z = 1 + j*b;
